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Abstract 

The predictions of an improved standard solar model are compared with the obser- 
vations of the four solar neutrino experiments. The improved model includes premain 
sequence evolution, element diffusion, partial ionization effects, and all the possible nu- 
clear reactions between the main elements. It uses updated values for the initial solar 
element abundances, the solar age, the solar luminosity, the nuclear reaction rates and 
the radiative opacities. Neither nuclear equilibrium, nor complete ionization are as- 
sumed. The calculated *B solar neutrino flux is consistent, within the theoretical and 
experimental uncertainties, with the solar neutrino flux measured by Kamiokande. The 
results from the ^^Cl and ^^Ga radiochemical experiments seem to suggest that the 
terrestrial ^Be solar neutrino flux is much smaller than that predicted. However, the 
present terrestrial "defecit" of ^Be solar neutrinos may be due to the use of inaccurate 
theoretical neutrino absorption cross sections near threshold for extracting solar neu- 
trino fluxes from production rates. Conclusive evidence for a real deficit of ^Be solar 
neutrinos will require experiments such as BOREXINO or HELLAZ. A real defecit of 
^Be solar neutrinos can be due to either astrophysical reasons or neutrino properties 
beyond the standard electroweak model. Only future neutrino experiments, such as 
SNO, Superkamiokande, BOREXINO and HELLAZ, will be able to provide conclusive 
evidence that the solar neutrino problem is a consequence of neutrino properties beyond 
the standard electroweak model. Earlier indications may be provided by long baseline 
neutrino oscillation experiments. 

1 Introduction 

The sun is a typical main sequence star which generates its energy by fusion of protons 
into Helium nuclei through the pp and CNO nuclear reactions chains which also produce 
neutrinos. These neutrinos have been detected on Earth in four pioneering solar neutrino 
experiments, the radiochemical Chlorine experiment at Homestake (Cleveland et al. 1995 
and references therein), the electronic light water Cerenkov experiment at Kamioka (Tot- 
suka, these proceedings) and the two radiochemical Gallium experiments, GALLEX at Gran 
Sasso (Kirsten, these proceedings) and SAGE at the Baksan (Gavrin, these proceedings). 
They provide the most direct evidence that the sun generates its energy via fusion of Hydro- 
gen into Helium. However, it has been claimed (e.g., Bahcall 1995, Hata et al. 1994) that all 



four experiments measured solar neutrino fluxes significantly smaller than those predicted 
by standard solar models (SSM) (e.g., Bahcall and Ulrich 1988; Turck Chieze ct al. 1988; 
Sackman et al. 1990, Bahcall and Pinsonneault 1992 (BP92), Turck-Chieze and Lopes 1993 
(TL93); Castellani et al. 1994; Kovetz and Shaviv 1994; Christensen - Dalsgaard 1994; Shi 
et al. 1994; Bahcall and Pinsonneault 1995 (BP95); See, however, Dar and Shaviv 1994 
(DS94); Shaviv (1995): Dzitko et al. 1995; Dar and Shaviv 1995 (DS95). 

This discrepancy has become known as the solar neutrino problem. Three types of 
solutions to the solar neutrino problem have been investigated: 

(a) Experimental Solution: Perhaps the accuracy of the results of the solar neutrino 
experiments has been overestimated and unknown systematic errors are largely responsible 
for the solar neutrino problem. 

(b) Astrophysical Solutions: Perhaps the standard solar models do not provide suffi- 
ciently accurate description of the present day sun and/or the neutrino producing reactions 
in the sun. 

(c) New Physics Solutions: Perhaps new electroweak physics beyond the standard elec- 
troweak model is responsible for the solar neutrino problem. 

The chances that possibility (a) is responsible for the solar neutrino problem have been 
greatly reduced by the GALLEX Chromium source experiment performed last year (Ansel- 
mann et al. 1995). This experiment is the first full demonstration of the reliability of the 
radiochemical technique for the detection of solar neutrinos. In particular, it excludes the 
possibility of any unidentified important sources of systematical errors, such as hot atom 
chemistry, in the radiochemical experiments. 

Astrophysical solutions have been ruled out by many authors on the grounds that the 
combined results from the four solar neutrino experiments show a large deficit of ^Be solar 
neutrinos relative to that expected from the ^B solar neutrino observations in Kamiokande 
(the ratio between proton capture in the sun by ^Bc, which produces ^B neutrinos, and 
electron capture by ^Be, which produces ^Be neutrinos, is almost model independent): 

1) The CI experiment with an energy threshold lower than that of Kamiokande measured 
a neutrino capture rate smaller than that expected from the ^B solar neutrino flux measured 
in Kamiokande, leaving no room for a significant ^Bc solar neutrino flux. 

(2) The Gallium experiments measured a neutrino capture rate which is consistent with 
the solar luminosity only if the sun does not produce a signiflcant flux of ^Be neutrinos. 

Consequently, it has been concluded (see e.g., Bahcall these proceedings) that non stan- 
dard neutrino properties, such as neutrino flavour mixing, are required in order to explain 
the solar neutrino observations. 

However, I do not consider yet the ^Be solar neutrino problem as a solid evidence for new 
electroweak physics, as claimed by various authors (e.g., Bahcall and Bethe 1993; Bludman 
et al. 1993; Castellani et al. 1994; Hata et al. 1994; Berezinsky 1994; Kwong and Rosen 
1994; Bahcall 1994; Parke 1995; Hata and Langackcr 1995). This is because the standard 
solar models are only approximate and simplified descriptions of the real and complex sun, 
because of very little experimental and theoretical knowledge of dense plasma effects on 
nuclear reaction rates and on energy transport at solar conditions, and because standard 
physics solutions are still possible. In particular, below I will argue that: 
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i) The solar neutrino flux observed by Kamiokande is consistent with the predictions 
of the standard solar models of Dar and Shaviv (1995) and of Dzitko et al. (1995). 

ii) Suppression of the neutrino capture rates near threshold in ^''Cl and in ^^Ga may 
be responsible for the discrepancy between the observed capture rates of solar neutrinos in 
^^Cl and in ^^Ga and those calculated from the standard solar model. 

2 The Standard Solar Model - An Overview 

The standard solar model (e.g. Bahcall 1989 and references therein) is a physical description 
of the sun based on the standard stellar evolution equations (e.g., Clayton 1968), which are 
used to calculate its evolution from its premain sequence Hayashi phase to its present state. 
The model assumes a complete spherical symmetry, no mass loss or mass accretion, no an- 
gular momentum gain or loss, no differential rotation and a zero magnetic field through the 
whole solar evolution . It uses best available input physics ( equations of state, nuclear cross 
sections, radiative opacities, condensed matter effects) and the following initial conditions: 

1. Fully convective, homogeneous, spherically symmetric protostar. 

2. Initial mass of AfQ = 1.99 x 10-^^ gm. 

3. Zero angular momentum, no differential rotation, no magnetic field. 

4. Initial chemical composition deduced from primitive meteorites, the solar photosphere, 
the solar wind, the lunar soil, the interstellar medium and the photospheres of nearby 
stars. 

The calculations are iterated, treating the unknown initial '^He abundance and the mix- 
ing length in the convective zone (roughly the size of the pressure scale height) as adjustable 
parameters, until the present observed properties of the sun are reproduced at its present 
estimated age, = 4.57 GY. These include (Particle Data Group 1994): 

1. The solar luminosity Lq = 3.844 x 10^^ erg ■ sec~^. 

2. The solar radius Rq = 6.9599 x 10^°cm. 

3. The observed solar surface element abundances. 

4. Internal structure consistent with helioseismology. 

The output of the calculations includes the present-day density profile p(r), temperature 
profile T(r) and chemical composition profile [Xj(r)] of the sun. They can be compared with 
information extracted from helioseismology (see Christensen-Dalsgaard, these proceedings). 
They can also be used to calculate the expected fluxes of solar neutrinos which are produced 
in the fusion of hydrogen into deuterium {pp De^Ve and pep DUf,), in electron capture 
by ^Be and in /? decay of ^B, ^^N, ^^O and ^^F using the standard electroweak theory and 
the best available nuclear reaction/decay rates. 
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Our calculations were performed with an updated version by Dar and Shaviv (1994; 
1995) of the solar evolution code of Kovetz and Shaviv (1994). This updated version includes 
premain sequence evolution, diffusion of all elements, partial ionization effects, and all the 
significant nuclear reactions between the various elements which the sun is made of. It uses 
updated values for the initial solar element abundances, the solar age, the solar luminosity, 
the nuclear reaction rates and the radiative opacities. Neither nuclear equilibrium, nor 
complete ionization are assumed. It employs a very fine zoning of the sun and accurate 
numerical procedures to integrate the solar evolution equations from zero age until the 
present day. We refer the reader to its detailed description by Kovetz and Shaviv (1994) 
and by Dar and Shaviv (1995). Here we highlight only the choice of initial solar composition 
and nuclear reaction rates. 

2.1 Initial Chemical Composition 

The initial element abundances influence significantly solar evolution and the present den- 
sity, chemical composition and temperature in the solar core, which determine the solar 
neutrino fluxes. In particular, the calculated radiative opacities, which determine the tem- 
perature gradient in the solar interior, are very sensitive to the abundances of the heavy 
elements (which are not fully ionized in the sun). 

Four major sources of information on the initial solar abundances were used: The chem- 
ical composition of the most primitive class of meteorites (type I carbonaceous chondrites), 
the solar photospheric abundances, the chemical composition of the solar wind and the 
lunar soil. 

We have assumed that the initial solar heavy metal abundances are given approximately 
by the meteoritic (carbonaceous CI chondrites) values of Grcvcssc and Noels (1993a). They 
better represent the initial values than the photospheric ones, which have changed by mixing, 
convection and diffusion during the solar evolution, and which are known with much less 
accuracy. 

The meteoritic CNO abundances are very much different from those found in the solar 
photosphere. In particular. Carbon and Nitrogen are underabundent in CI chondrites by 
about an order of magnitude compared with their photospheric abundances (e.g., Sturen- 
berg and Howlweger 1990; Grevesse et al. 1990). Therefore, we have used initial CNO 
abundances which yield their present photospheric values when diffusion is included in the 
solar evolution calculation. 

The initial ''He mass fraction in the solar nebula which is known only approximately, 
0.24 <Y< 0.30, has been treated as an adjustable parameter. 

The initial solar abundances of ^He and D, were taken to be pHe]/[H]=(1.5±0.3) x 10~^ 
and [D]/[H] = (2.6 ± 1.0) x 10"^ respectively. These values were inferred by Gciss (1993) 
from analysis of solar wind particles captured in foils exposed on the moon, from the lunar 
surface and from studies of primitive meteorites. 

The photospheric abundances of ^Li, ^Bc and ^^B are smaller by a factor of nearly 150, 3 
and 10, respectively, than their meteoritic abundances. The origin of such large differences 
is not clear. However, the initial solar (meteoritic) abundances of these elements are very 
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small and do not play any significant role in solar evolution (their solar depletion, however, 
may store important information on the true history of the sun). 

2.2 Nuclear Reaction Rates 

The uncertainties in the nuclear reaction rates at solar conditions arc still large due to: 
(1) uncertainties in the measured cross sections at laboratory energies, (2) uncertainties in 
their extrapolations from laboratory energies down to solar energies, (3) uncertainties in 
dense plasma effects (screening, correlations and fluctuations) on reaction rates. Neither the 
microscopic methods (for reviews see, e.g., Langanke 1991) such as the Resonating Group 
Method (RGM) and the Generator Coordinate Method (GCM), nor the potential models 
such as the Optical Model (OM) and the Distorted Wave Born Approximation (DWBA), 
can predict accurately enough the low energy cross sections. Consider for instance the 
reaction ''Be(p,7)^B. The RGM and the GCM, which are currently considered to be the best 
theorerical methods for calculating direct nuclear reactions, predict (see, e.g. Descouvemont 
and Baye 1994, Johnson et al. 1992) 517(0) pa 25 — 30 • b. However, a simple inspection 
of their predictions reveals that they poorly reproduce the magnitude of the measured 
cross section, the position of the resonance, the width of the resonance, the height of the 
resonance and the observed shape of the cross section as function of energy. To avoid these 
discrepancies only the energy dependence of these models has been used to extrapolate 
the measured cross sections to E = 0, yielding (see, e.g., Johnson et al. 1992) £'17(0) ~ 
22.4 eV - b . This value has been used in BP92 and BP95. However, the procedure used by 
Johnson et al. (1992) is rather ad hoc and questionable: (a) Their model does not reproduce 
the measured energy dependence of the cross section at lab energies, (b) The "average 
cross section" which they extrapolated to solar energies was obtained by averaging cross 
sections which differ by many standard deviations and have different energy shapes (the 
cross sections measured by Kavanagh (1960) and by Parker (1968) differ by more than 3a 
from the cross sections measured later by Vaughn (1970) and by Filippone (1983) in the same 
energy ranges; [Kavanagh] /[Filippone] = 1.34 it 0.11, [Parker] /[Vaughn] = 1.42 ±0.13; see, 
e.g., Gai 1995). 

Dar and Shaviv (1994) noted that sub-Coulomb radiative captures and transfer reactions 
take place mainly when the colliding nuclei are far apart. They argued that since optical 
models describe well the shapes of the bound state and relative motion wave functions 
outside the nuclear potential, they should be preferred for extrapolating the laboratory 
cross sections to solar energies. Alternatively, they proposed to remove the trivial energy 
dependence due to Coulomb barrier penetration to the effective distance R (calculated 
or best fitted) where the reaction takes place and then to use a simple polynomial fit to 
extrapolate the measured cross sections to solar energies. The resulting astrophysical S 
factors, which were obtained this way, are summarized in Table II. They include: 

817(0): Extrapolations of the cross sections measured by Vaughn (1970) and by Fil- 
ippone(1983), using either simple potential models or the very general properties of sub- 
Coulomb cross section, gave 6*17(0) ~ 17 eV ■ b (Barker and Spear 1986; Dar and Shaviv 
1994; Kim et al. 1994). Similar values were obtained also by different types of experiments: 
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Analysis of the virtual photodissociation reaction 7t,+^B^ p+^Bc in the Coulomb field of 
Lead (Motobayashi et al. 1994) gave S'i7(0) f« 17 eV -b. A similar value, S'i7(0) ^ 17.6 eV -b, 
was estimated from the virtual reaction py +^ Be ^ B measured through the proton trans- 
fer reaction ^He+^Be— >D-|-^B (Xu et al. 1994). Consequently, we a have adopted the value 
Sir = 17 eV-b in our standard solar model calculations rather than the value Sn = 22.4 eV-b 
which was used in BP92 and BP95. 

834(0): The value >S'34(0) = 0.51 ± 0.02 keV ■ b was obtained from measurements of 
the prompt 7-ray emitted in the reaction ^He +^ He -^"^ Be + 7, while measurements of 
the induced ^Be activity led to a weighted average 534(0) = 0.58 ± 0.02 keV ■ b which 
is different by 3.5 standard deviations. The origin of this discrepancy is still not known 
(Hilgemeier et al. 1988). Normalization to known cross sections favour the lower value. 
Consequently, the value 534(0) = 0.45 keV ■ b was obtained by using the energy depen- 
dence of the data of Krawinkel et al., 1982 to extrapolate the measured cross sections for 
^He +^ He -^'^ Be + 7 from prompt gamma ray emission (Parker and Kavanagh 1963; Os- 
borne et al. 1982; Krawinkel et al. 1982 (multiplied by 1.4); Hilgemeier et al. 1988) to solar 
energies. Note that the 534(0) = 0.533 keV ■ b was used in BP92 and ^34(0) = 0.524 keV ■ b 
was used in BP95. 

833(0): Extrapolation of the low energy data of Greife et al. (1994), Krauss et al. 
(1987) and Dawarakanath and Winkler (1971) on the reaction ^He -|-^ He -^^ He -|- 2p gave 
533(0) = 5.6 MeV ■ b. Essentialy the same value, 533(0) = 5.57 MeV ■ b, was obtained by 
Krauss et al. (1987) and by Greife et al. (1944) by applying a polynomial fit to their data. 
The value 533(0) = 5.0 MeV ■ b was used in BP92 and 533(0) = 4.99 MeV ■ b was used in 
BP95. 

Sii(O): The reaction p + p^D + e+ + fe has a cross section which is too small to be 
measured directly in the laboratory. The weak isospin related reactions Pe+D— > e^-l-n-l-n, 
z^e+D^ Pg+P+n; and n + p — > Z? + 7 have been measured and can be used to obtain the 
relevant nuclear matrix element needed for calculating the cross section for p+p^D+e~^ +i'e- 
This procedure yields a best value Sn (0) ~ 4.07 x 10~^^ keV ■ b which is consistent with 
the value used by Caughlan and Fowler (1988). It is 4.2% larger than the value 5ii(0) ^ 
3.896 X 10^^^ keV ■ b calculated recently by Kamionkowski and Bahcall (1994). 

Screening Enhancement of Reaction Rates: Screening of target nuclei by electrons 
known to enhance significantly laboratory nuclear cross sections at very low energies (e.g., 
Engstler et al. 1988). A complete theoretical understanding of the effect is still lacking (see, 
e.g., Shoppa et al. 1993; Rolfs 1994 and references therein). The screening enhancement 
factors of the nuclear reaction rates near the center of the sun (Bahcall 1989) are quite 
considerable, being 5%, for the pp and pep reactions, 20% for the ^He^He, ^He'^He and 
p^Be reactions, and 30%, 35% and 40% for the p capture by the C, N and O isotopes, 
respectively. However, these screening enhancement factors are based on the Debye-Huckel 
approximation which is not valid for the conditions prevailing near the center of the sun. 

Reliable estimates of the nuclear reaction rates in dense plasma require numerical N 
body simulations of the behaviour of electrons and ions in a dense plasma. Such simula- 
tions can be used to evaluate the effects of screening, correlations and fluctuations on the 
thermonuclear reaction rates in dense plasmas (Shaviv and Shaviv 1995). They have not 
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yet been incorporated in our solar evolution code. 

However, to test the sensitivity of the standard solar models to the screening corrections 
we have carried out the calculations with, and without, the standard screening enhancement 
factors of all the thermonuclear reaction rates. We have found (Dar and Shaviv 1994) that 
removing/including the screening enhancement factors for all nuclear reaction rates had only 
a small net effect on the calculated solar neutrino fluxes, due to accidental cancellations. 
Screening factors may, however, play important role if their ratios for the different reactions 
are changed. 

2.3 Radiative Opacities 

The radiative opacities depend on the local chemical composition, density and temperature 
in the sun. We have used radiative opacity tables computed by the OPAL group at Lawrence 
Livermore National Laboratory that were kindly provided to us by F.J. Rogers. They are 
updated version of the OPAL tables of Rogers and Iglesias (1992) for the most recent 
determination by Grevesse and Noels (1993a) of the heavy element composition of the sun 
from the meteoritic and photospheric data. 

2.4 Diffusion 

The Kovetz-Shaviv code calculates the diffusion of all the individual elements from the 
premain sequence phase to the present age. The binary and thermal diffusion coefficients 
depend on the squared ionic charges. The ionization state of each element in every shell 
is calculted by solving the Saha equations for all the elements in each shell. All elements 
with mass fractions less than 10~^ are treated as trace elements. Their collisions with other 
trace elements are neglected. 

As a consequence of diffusion the surface and internal element abundances change for 
each element in a slightly different way. Diffusion depletes the surface abundances of ^He 
and the heavy elements. For C, N and O only their present photospheric abundances arc 
known. Consequently, in our solar model calculations their initial solar abundances were 
adjusted to reproduce their observed photospheric abundances. For all other elements the 
initial meteoritic abundances were used to predict their final surface abundances which can 
be compared with their observed photospheric abundances. The unknown initial abundance 
of ^He was treated as a free parameter. It was adjusted to best reproduce the presently 
observed sun. Its predicted surface abundance today can be compared with the value 
derived from helioseismology (Hernandez and Christensen - Dalsgaard 1994). 

3 Results and Comparison With Observations 

Our predictions of the solar neutrino fluxes for three standard solar models are summarized 
in Table Ilia and are compared with the results from the four solar neutrino experiments. 
The three models differ only in their teatment of element diffusion: The model labeled 
DSND does not include element diffusion. The models labeled DS94 and DS95 include 
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element diffusion. Model DS94 assumes that the initial heavy metal abundances in the sun 
were equal to their meteoritic values, and those of C, N and O were equal to their observed 
photospheric abundances, as summarized in Table I. In model DS95 the initial abundances 
of C, N and O were adjusted to yield their present day photospheric abundances while 
the heavy metal abundances were assumed to be those found in primitive meteorites, as 
summarized in Table I. Table Illb presents some physical characteristics of the three models. 
As can be seen from Table Ilia, all three models yield a solar neutrino flux consistent 
with that measured by Kamiokande. However all three models predict capture rates in the 
Chlorine and Gallium experiments which are significantly larger than those measured by 
Homestake, GALLEX and SAGE (we have used the neutrino cross sections from Table 8.2 
of Bahcall (1989) to convert solar neutrino fluxes to capture rates). 

Table IV present a comparison between four solar models. The model labeled BP92 is 
the best model of BP92. It includes diffusion of protons and ^He but not of other elements. 
The model labeled BP95 is the best model of BP95 which includes also diffusion of the heavy 
metals but assumes that all the heavy elements diffuse like fully ionized iron. The predictions 
of the DS models differ significantly from the BP models because of different input physics, 
approximations and numerical methods. Most of the differences are due to the use of 
different reaction rates as summarized in Table II and the different treatments of diffusion. 
This is demonstrated in Table V where we present a comparison between the best model of 
BP95 without diffusion labeled BP95ND, and a solar model, labeled DS(BPND) calculated 
with the Kovetz-Shaviv stellar evolution code with the same physical and astrophysical input 
parameters and the same nuclear reaction rates used in BP95ND, and without inclusion of 
element diffusion. As can be seen from Table V the two calculations yield similar results. 
Even the fluxes of ^B and CNO solar neutrinos which, under the imposed solar boundary 
conditions, are very sensitive to the central solar temperature differ by less than 4%. These 
remaining differences are probably due to the use of different equations of state, numerical 
methods, fine zoning and time steps in the two codes and due to the inclusion of premain 
sequence evolution in our code. 

To emphasize the important role that might be played by diffusion. Table V also include 
the current best solar models of Dar and Shaviv (1995) and of Bahcall and Pinsonneault 
(1995), which include diffusion of all elements. As can be seen from Tables V, Bahcall and 
Pinsonneault (1995) found rather large increases in their predicted ^Be, ^B, ^^N, ^^O and 
^^F solar neutrino fluxes; 14%, 36%, 52%, 58%, and 61%, respectively, compared with their 
model (BP95ND) with no diffusion. These induce 36%, 33%, 9% increases in the predicted 
rates in Kamiokande, Homestake, and in GALLEX and SAGE, respectively. However, we 
predict more moderate increases due to diffusion, 4%, 10%, 23%, 24% and 25%, respectively, 
in the above fluxes, which correspond to 10%, 10% and 2% increases in the predicted rates in 
Kamiokande, Homestake, and in GALLEX and SAGE, respectively. The differences in the 
effects of diffusion in DS94 and BP95 are mainly due to two reasons: (a) In the calculations 
of Bahcall and Pinsonneault (BP95) all heavy elements were assumed to diffuse at the same 
rate as fully ionized iron while in the Dar-Shaviv calculations (DS94) followed the diffusion 
of all the elements separately and used diffusion coefficients for the actual ionization state 
of each element, (b) Bahcall and Pinsonneault assumed that the meteoritic abundances 
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represent the solar surface abundances today and not their initial values. They adjusted 
their initial values to reproduce surface abundances today equal to the meteoritic values. 
Thus they have actually used an initial ratio Z/X = 0.0285 (see Table IV of BP95) while 
the observed ratio in meteorites is ZjX = 0.0245 (Grevesse and Noels 1993). Dar and 
Shaviv (1994) used the meteoritic values for the initial metalic abundances and predicted 
present day depleted surface abundances. Unfortunately, the uncertainties in the measured 
photospheric CNO abundances (typically 30%-40%) are much larger than their predicted 
depletion (typically 10%) and do not allow a reliable test of the predicted photospheric 
abundances. 

Standard solar models do not explain the large depletion in the abundance of Li in the 
solar photosphere compared to its meteoritic abundance. Such a depletion may indicate 
significant mixing of surface material to solar depths where Li is burned efficiently. Mixing 
may , however, inhibit the inward diffuson of the heavy elements. Thus, it is not clear 
whether solar models with diffusion provide a more realistic description of the sun than 
solar models without diffusion, although helioseismology data is better explained by SSM's 
with diffusion. See, e.g., Cristensen-Dalsgaard et al. 1993. 

4 Where Are The ^Be Neutrinos ? 

Electron capture by ^Be into the ground state of ^Li produces 862 keV neutrinos. The 
threshold energy for neutrino absorption by ^^Cl is 814 keV. Thus, absorption of ^Be neu- 
trinos by ^^Cl produces 48 KeV electrons. The average energy of the pp solar neutrinos is 
265 KeV. The threshold energy for neutrino absorption in '''^Ga is 233 KeV. Consequently, 
the produced electron has a typical energy of 33 keV. The de Broglie wave lengths of such 
electrons are larger than the Bohr radii of the atomic K shells in CI and Ga and their ener- 
gies are similar to the kinetic energies of electrons in the K shells. Consequently, screening 
of the nuclear charge by atomic electrons and final state interactions (exchange effects, ra- 
diative corrections, nuclear recoil against the electronic cloud, etc.) may slightly reduce the 
absorpton cross sections of pp neutrinos in ^^Ga (perhaps making room for the expected 
contribution of ^Be in Gallium ?) and of ^Be neutrinos in ^^Cl (perhaps making the solar 
neutrino observations of Kamiokande and the Homestake experiment compatible). It is 
interesting to note that although final state interactions in Tritium beta decay have been 
studied extensively, they do not explain yet why its end-point spectrum [Ef. ~ 18.6 kEV) 
yields consistently, in all recent measurements, a negative value for the squared mass of the 
electron neutrino. Final state interactions in ^^Cl and '''^Ga are expected to be much larger 
because of their much larger values of Z. Note also that the above explanation implies that 
experiments like BOREXINO and HELLAZ will observe the full ''Be solar neutrino flux 
while the MSW solution predicts that it will be strongly suppressed. 

Even if the ^Be solar neutrino flux is strongly suppressed, it does not elliminate yet 
standard physics solutions to the solar neutrino problem. For instance, collective plasma 
physics effects, such as very strong magnetic or electric fields near the center of the sun, 
may polarize the plasma electrons, and affect the branching ratios of electron capture by 
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^Be (spin 3/2") into the ground state (spin 3/2", E^,^ = 0.863 MeV, BR=90% and the 
excited state (spin 1/2^, Ey^^ = 0.381 MeV, BR=10%) of ^Li. Since solar neutrinos with 
-Ei/g = 0.381 MeV are below the threshold (0.81 MeV) for capture in ^'^Cl and have a capture 
cross section in ^^Ga that is smaller by about a factor of 6 relative to solar neutrinos with 
E,^^ = 0.863 MeV, therefore a large suppression in the branching ratio to the ground state 
can produce large suppressions of the ^Be solar neutrino signals in ^^Cl and in '^^Ga. 

5 Conclusions: 

Solutions to the solar neutrino problem which do not invoke physics beyond the standard 
electroweak model are not ruled out yet: 

The solar neutrino problem may be a terrestrial problem. The neutrino capture cross sec- 
tions near threshold in the radiochemical experiments may be different from the calculated 
cross sections. The inferred solar neutrino fluxes from the GALLEX and HOMESTAKE 
experiments may be different from the true solar neutrino fluxes. They do not established 
beyond doubt that the ^Be solar neutrino flux is strongly suppressed. BOREXINO and 
HELLAZ may be able to do that. 

The solar neutrino problem may be a problem of the standard solar model. Namely, 
the model may not provide yet an accurate enough description of the sun and the nuclear 
reactions that take place in it. The deviations of the experimental results from those 
predicted by the standard solar models may reflect the approximate nature of the standard 
solar models (which neglect, or treat only approximately, many effects and do not explain 
yet solar activity nor the surface depletion of Lithium, Berilium and Boron relative to 
their meteoritic values, that may or may not be relevant to the solar neutrino problem). 
Improvements of the standard solar model should continue. In particular, dense plasma 
effects on nuclear reaction rates and radiative opacities, which are not well understood, 
may strongly affect the SSM predictions and should be further studied, both theoretically 
and experimentally. Relevant information may be obtained from studies of thermonuclear 
plasmas in inertial confinement experiments. Useful information may also be obtained from 
improved data on screening effects in low energy nuclear cross sections of ions, atomic beams 
and molecular beams incident on a variety of gas, solid and plasma targets. 

Better knowledge of low energy nuclear cross sections is badly needed. Measurement 
of crucial low energy nuclear cross sections by new methods, such as measurements of the 
cross sections for the radiative captures p +^ Be — >^ B + 7 and ^He +^ He — Be + 7 by 
photodissociation of ^B and ^Be in the coulomb field of heavy nuclei are highly desireable. 

The ^"^Ar production rate in ^^Cl may indeed be smaller than that expected from the 
total solar neutrino flux measured by electron scattering in the Kamiokande experiment. In 
that case neutrino oscillations, and in particular the MSW effect, may be the correct solution 
to the solar neutrino problem. Only future experiments like SNO and Superkamiokande will 
be able to supply a definite proof that Nature has made use of this beautiful effect. 
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Table I: Summary of Information on Abundances of Various Elements relative to Hy- 
drogen ( A=log([A]/[H])+12) in Primitive Meteorites, in the Solar Photosphere, in the Solar 
Wind and in the Local Interstellar Medium, Used In The DS Standard Solar Models. 



Element 


Abundance 


Source 


Reference 


D 


7.22 ± 0.05 


LISM, Meteorites 


Linsky 1993, Geiss 1993 


''He 


7.18 ± 0.08 


Meteorites, Solar Wmd 


Geiss 1993 


'Li 


1.54 ± O.OX 


Meteorites 


Anders and Grevesse 1989 


^Be 


1.13 ± O.OX 


Meteorites 


Anders and Grevesse 1989 




8.55 ± 0.05 


Photosphere 


Grevesse and Noels 1993 




6.60 ± 0.05 


Photosphere 


Grevesse and Noels 1993 




7.97 ± 0.07 


Photosphere 


Grevesse and Noels 1993 


ibQ 


8.78 ± 0.07 


Photosphere 


Grevesse and Noels 1993 


"'"Ne 


8.08 ± 0.06 


Photosphere 


Grevesse and Noels 1993 




6.33 ± 0.03 


Meteorites and Photosphere 


Grevesse and Noels 1993 




7.58 ± 0.05 






97 \ 1 


6.47 ± 0.07 






^'^Si 


7.66 ± 0.05 






31p 


5.45 ± 0.04 








7.21 it 0.06 






•^^Cl 


5.5 ± 0.3 






Ar 


0.52 ± 0.1 






40Ca 


6.36 ±0.02 






40k 


4.85 






45Sc 


3.08 






48Ti 


5.02 ± 0.06 






SOy 


3.99 






52Cr 


5.67 ±0.03 






55Mn 


5.39 ± 0.03 






56pg 


7.50 ± 0.04 






^^Cu 


4.15 






58Ni 


6.25 ±0.04 






64Zn 


4.33 






z/x 


0.0245 
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Table II: Comparison Between the Astrophysical S Factors for the pp-chain Reactions 
used in BP95 and in DS94 and DS95. The values of S are given in keV ■ b Units. 



Reaction 


SBP(O) 


S^^(O) 


iH(p,e+i/e)D 


3.896 X 10"^^ 


4.07 X 10 


iH(pe-z/e)D 


Bahcall 89 


CF88 


3He(3He, 2p)He^ 


4.99 X 10^ 


5.6 X 10^ 


3He(^He,7)'^Be 


0.524 


0.45 




Bahcall 89 


CF88 


^Be(p,7)'B 


0.0224 


0.017 



Table Ilia: Comparison Between Solar Neutrino Fluxes Predicted by the current best 
Standard Solar Models of Dar and Shaviv, with and without element diffusion, and the 
Solar Neutrino Observations. 



v Flux 






DSND 


DS94 


DS95 


Observation 


Experiment 


<pu{pp) [10^'' cm~ 






6.10 


6.06 


6.03 






(puipap) [lO^cm" 


-2.,- 




1.43 


1.42 


1.40 












4.03 


4.00 


4.20 


« (/.ssM(7Be) 


ALL 




^s- 




2.54 


2.60 


2.87 


2.9 ±0.4 


Kamiokande 




-2 s 




3.21 


3.30 


3.94 






-2s- 




3.13 


3.19 


3.88 






M^'^F) [lO^cm- 


-2s- 




3.77 


3.84 


4.71 






^{(l>a)ci [SNU] 






4.2 ± 1.2 


4.2 ± 1.2 


4.6 ± 1.6 


2.55 ± 0.25 


Homestake 


E(0C7)Ga [SNU] 






116 ±6 


116 ±6 


119 ±7 


79± 12 


GALLEX 


S(<^a)Ga [SNU] 






116 ±6 


116 ±6 


119 ±7 


74± 16 


SAGE 



16 



Table Illb: Characteristics of the DS Solar Models in Table Ilia (c=center; s=surface; 
bc=base of convective zone; N = log{[N]/[H]) + 12). 



Parameter 


DSND 


DS94 


DS95 




1.553 


1.554 


1.562 


Pc [a 


154.9 


155.3 


155.9 




0.3491 


0.3462 


0.3381 


Yc 


0.6333 


0.6359 


0.6425 




0.01757 


0.01802 


0.01940 




0.6978 


0.7243 


0.7171 




0.2850 


0.2597 


0.2658 


Zs 


0.01703 


0.01574 


0.01696 




8.55 


8.50 


8.55 




7.97 


7.92 


7.97 




8.87 


8.82 


8.87 




8.08 


8.03 


8.08 


Xs{>^^ Mg) 


0.00464 


0.00414 


0.00415 


Rconv [R/Rq] 


0.7306 


0.7105 


0.7066 


The [lO^K] 


1.97 


2.10 


2.12 


Teff [K] 


5895 


5920 


5919 



Table IVa: Comparison Between Solar Neutrino Fluxes Predicted by the Dar-Shaviv 
Models and the Bahcall-Pinsonneault best Solar Models. 







BP92 


(t)u{vp) [lO^'^ cm 




6.00 


(f)u{pep) [10^cm~ 


-2.,-!] 


1.43 






4.89 






5.69 


M^^N) [lO^cm 


-2^-1] 


4.92 






4.26 


M^'^F) [lO^cm" 




5.39 


^{(Pa)ci [SNU] 




8±1 


^{ct>a)Ga [SNU] 




132 ±7 



DS94 


BP95 


DS95 


6.06 


5.91 


6.03 


1.42 


1.40 


1.40 


4.00 


5.15 


4.20 


2.60 


6.62 


2.87 


3.30 


6.18 


3.94 


3.19 


5.45 


3.88 


3.84 


6.48 


4.71 



4.2 ±1.2 9.3 ±1.4 4.6 ±1.6 
116 ±6 137 ±8 119 ±7 
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Table IVb Characteristics of the BP95, DS94, and DS95 Solar Models in Table Ilia 
(c=center; s=surface; bc=base of convective zone; N = log([N]/[H]) + 12). 



Parameter 


BP95 


DS94 


DS95 


Tc [W'^K] 


1.584 


1.554 


1.562 


Pc [g cm 


156.2 


155.3 


155.9 




0.3333 


0.3462 


0.3381 


Yc 


0.6456 


0.6359 


0.6425 




n no 1 1 

0.0211 


0.01802 


0.01940 




0.7351 


0.7243 


0.7171 




0.2470 


0.2597 


0.2658 




0.01798 


0.01574 


0.01696 




8.55 


8.50 


8.55 




7.97 


7.92 


7.97 




8.87 


8.82 


8.87 




8.08 


8.03 


8.08 


X,(>24 Mg) 




0.00414 


0.00415 


Rconv [R/Rq] 


0.712 


0.7105 


0.7066 


Tbc [lO^K] 


2.20 


2.10 


2.12 


Teff [K] 




5920 


5919 



Table V: Comparison between the solar neutrino fluxes calculated from the best stan- 
dard solar model with no diff'usion of BP95 and those calculated with the Dar-Shaviv SSM 
code with the same nuclear reaction rates, opacities, composition and astrophysical param- 
eters. The predictions of the current best standard solar models of Dar and Shaviv and of 
Bahcall and Pinsonneault are also included. 







BP95(ND) 


DS(BPND) 


BP95 


DS95 


(f>u{pp) [lO-'^^cm 




6.01 


6.08 


5.91 


6.03 


(t)y{pep) [10*cm~^s" 




1.44 


1.43 


1.40 


1.40 


(t)uCBe) [lO^cm-^s 




4.53 


4.79 


5.15 


4.20 


(i)S{^B) [lO^cm-^s- 




4.85 


5.07 


6.62 


2.87 


(l)S{^^N) [lO^cm-^s 
4)S{^^0) [lO^cm-^s- 




4.07 


2.50 


6.18 


3.94 




3.45 


3.38 


5.45 


3.88 


(t^ul^^F) [lO^cm-^s- 




4.02 


4.06 


6.48 


4.71 


T.{(t)a)ci [SNU] 




7± 1 


7± 1 


9.3 ±1.4 


4.6 ± 1.6 


S(0a)Ga [SNU] 




127 ±6 


128 ±7 


137 ±8 


119 ±7 
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